UICK RESPONSE to relieve the disaster by supplying proper and critical material attracts a lot of research efforts in the world for its profound impact on social life and stability

Timely supply of vital materials to disaster hit areas plays a critical role in emergency relief. The problem involves warehouse selection, and fleet routing and scheduling so as to meet demand in the strict time window. The problem is NP-hard in general and extremely difficult to solve. The congestion caused by heavy traffic further aggravates the problem. To obtain a scalable solution, a new method based on successive subproblem solving in Lagrangian relaxation (LR) framework is developed. The route capacity and location selection constraints are relaxed by Lagrange multipliers, and the problem is converted into a two-level optimization problem. The subproblems at the lower level are solved successively in dual iterations with convergence assurance so that the indecomposable location constraints can be incorporated. A systematic method is developed to obtain a feasible solution by adding the once relaxed constraints back into the dual problem successively in feasibility iterations. Convergence proof of the new method and its properties are presented. Numerical results show that the new method is effective and efficient, and can be applied to large scale problems. Note to practitioners—The motivation of our paper is to simultaneously resolve the warehouse selection/location and the fleet routing and scheduling problem for emergency material delivery. Currently, most research work reported in the literature for emergency material supply problem focuses on either location selection problem or fleet routing and scheduling problem. In many practical applications, these two types of problems have to be synergistically considered and solved for timely material delivery. Since the problem is NP-hard, it is extremely difficult to solve by general integer programming methods. A new method based on the successive subproblem solving scheme is developed to solve this problem. The near optimal solution is obtained systematically by adding the once relaxed constraints back into the dual problem successively in feasibility iterations with convergence proof. The numerical testing results show that the new method can The research presented in this paper is supported in part by the National Natural Science Foundation (60921003, 60736027), 863 High Tech Development Plan (2007AA04Z154) and 111 International Collaboration Program, of China. Y. Han, X. Guan and Leyuan Shi are with Center for Intelligent and Networked Systems, Department of Automation, NLIST, Tsinghua University, Beijing, 100084 China; email: [email protected]; X. Guan is also with SKLMS Lab and MOE KLINNS Lab, Xian Jiaotong University, Xi’an, 710049 China; L. Shi is also with the Department of Industrial Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA. solve large scale problems with satisfactory duality gaps and is much more efficient than general integer programming methods.